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Abstract. The symmetric periodic Anderson model is well known to capture the essential physics of Kondo 
insulator materials. Within the framework of dynamical mean-field theory, we develop a local moment 
approach to its single-particle dynamics in the paramagnetic phase. The approach is intrinsically non- 
perturbative, encompasses all energy scales and interaction strengths, and satisfies the low-energy dictates 
of Fermi liquid theory. It captures in particular the strong coupling behaviour and exponentially small 
quasiparticle scales characteristic of the Kondo lattice regime, as well as simple perturbative behaviour in 
weak coupling. Particular emphasis is naturally given to strong coupling dynamics, where the resultant 
clean separation of energy scales enables the scaling behaviour of single-particle spectra to be obtained. 



PACS. 71.27.+a Strongly correlated electron systems; heavy fermions 
pounds and alloys; Kondo effect, valence fluctuations, heavy fermions 
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1 Introduction 

In the field of heavy fermion physics, the periodic An- 
derson model (PAM) has long played a key role in un- 
derstanding the rich range of behaviour arising in these 
lanthanide or actinide based compounds [1,2]. A natu- 
ral lattice generalization of the single-impurity Anderson 
model [2], each lattice site in the PAM contains a non- 
degenerate, correlated f-orbital that hybridises locally to 
a non-interacting conduction band. The model is certainly 
a simplification of material reality, albeit rather a good 
one for Ce-based systems where crystal field splittings re- 
duce the multiplicity of the f ^configuration to a Kramers 
doublet. Yet its simplicity is nominal: a microscopic un- 
derstanding of the PAM is far from complete, and even in 
low-temperature Fermi liquid phases many basic issues re- 
lating to the formation, nature and description of coherent 
quasiparticles remain open. 

In recent years, considerable progress in understanding 
correlated electron systems has been made within the pow- 
erful framework of dynamical mean-field theory (DMFT, 
reviewed in [3-6]); which is formally exact in the limit of 
infinite spatial dimensions, and provides a tangible local 
approximation in finite dimensions without trivialising the 
central role of interactions. Within DMFT the dynamics 
of the system become essentially local, and all correlated 
lattice-fcrmion models reduce to an effective quantum im- 
purity hybridizing self-consistently with the surrounding 
fcrmionic bath [3-6] . But therein lies an unsurprising dif- 
ficulty, for to solve such problems entails the ability to de- 
scribe an Anderson impurity model (AIM) with essentially 



arbitrary dynamics (w-dependence) in the hybridization 
function A e s(u>) that encodes the coupling between the 
fiducial impurity and the underlying host /bath. Such a 
description should be capable of (i) handling the problem 
on all energy scales and (ii) across the full range of (f- 
electron) interaction strengths U; and (iii) at low-energies 
must satisfy the dictates of Fermi liquid theory arising 
from adiabatic continuity to the non-interacting limit (as- 
suming such a state arises). As such it must be intrinsi- 
cally non-perturbative in order to capture the exponen- 
tially small scales symptomatic of strongly correlated be- 
haviour, yet must also recover perturbative, weak coupling 
behaviour and the non-interacting limit of the problem. 

The theoretical difficulties here are well known to be 
considerable, even for a pure AIM and naturally the more 
so for self-consistent lattice models; particularly in the 
context of dynamical properties such as single-particle ex- 
citation spectra, optical conductivities and associated trans- 
port properties. The PAM itself has of course been studied 
extensively within DMFT [3-6]: via numerical methods 
such as the numerical renormalization group (NRG) [7,8], 
quantum Monte Carlo (QMC) [9-12] and exact diagonal- 
ization (ED) [13]; and theoretical approaches that include 
perturbation theory (PT) in the interaction strength [14, 
15], iterated perturbation theory (IPT) [16,17], the lat- 
tice non-crossing approximation (LNCA) [18,19] and the 
simpler but related average t-matrix approximation [20], 
large- N mean-field theory [21,22], and the Gutzwiller ap- 
proach [23,24]. While much has been learned from these 
techniques, most suffer from well recognised limitations 
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and only the NRG [7,8] meets the above criteria. For ex- 
ample, finite-size effects render QMC and ED of limited 
value, the former being further confined to modest inter- 
action strengths and finite temperature (T). On the the- 
oretical side, perturbative methods are restricted to weak 
coupling and cannot recover exponentially small scales; 
NCA-based approaches violate Fermi liquid behaviour at 
low-energies and do not recover the non-interacting limit; 
while large-N mean-field theory, which cannot handle the 
full range of interaction strengths, amounts simply to a 
renormalization of the non-interacting limit and as such 
is applicable only on the very lowest-energy scales. 

The above comments underscore the need for the de- 
velopment of new, necessarily approximate theories. In 
the context of pure quantum impurity models (AIMs) we 
have recently initiated one such: the local moment ap- 
proach (LMA) [25-33], whose primary focus is dynamics 
and transport properties. The LMA satisfies the desider- 
ata listed above. It handles all energy scales on an equal 
footing, and its intrinsically non-perturbative nature en- 
ables it to capture the spin-fluctuation physics character- 
istic of the strong coupling Kondo regime, embodied dy- 
namically in the exponentially narrow, low-energy Kondo 
resonance; yet it also spans the full range of interactions, 
reducing asymptotically in weak coupling to straight sec- 
ond order PT in U [25,30,31]. Fermi liquid behaviour is 
moreover recovered at low energies - where appropriate, 
the latter point emphasising that the approach is not con- 
fined to the conventional Fermi liquid physics inherent to 
the metallic AIM, but can also handle the non-Fermi liq- 
uid behaviour and associated quantum phase transition 
arising e.g. in the pseudogap AIM [31-33]. Symmetric [25- 
29,31-33] and asymmetric [30] AIMs can now be handled 
within the LMA, which has also been extended to incorpo- 
rate finite-T [29] and the role of an applied magnetic field 
H [27,28]. Results for dynamics arising therefrom have 
also been shown [26,28-32] to give very good agreement 
with NRG calculations; and, for static magnetic proper- 
ties of the metallic AIM, with exact results provided by 
the Bethe ansatz [27,28]. 

In this paper we develop the LMA within DMFT to 
encompass the symmetric PAM [7-10,13,14,16,19], well 
known to be important in understanding the class of mainly 
cubic Kondo insulating materials such as CeNiSn, SmB 6 , 
Cc3Bi4Pt3 and YbBi2 (see e.g. [34,35]), whose low-energy 
electronic structure is characterised by an interaction-ren- 
ormalized hybridization gap. Specifically we consider here 
T = single-particle dynamics of the paramagnetic phase. 
After the relevant background (§2), the LMA for the lat- 
tice model is specified (§3), centering on the two-self-energy 
description and notion of symmetry restoration (§3.1) that 
underlie the approach. That discussion is general, appli- 
cable to an arbitrary diagrammatic approximation for the 
associated dynamical self-energies, E„{lo); the particular 
non-perturbative class of diagrams we implement here in 
practice is specified in §3.2. Following a brief discussion 
highlighting deficiencies of the static mean-field approxi- 
mation (§4), LMA results for single-particle dynamics are 
given in §5. As for the pure AIMs considered hitherto [25- 



33] , the LMA passes the criteria outlined above. Our pri- 
mary emphasis, albeit not exclusive, is naturally on the 
strong coupling regime of the PAM. Here, granted the abil- 
ity to capture exponentially small scales characteristic of 
the Kondo lattice, the resultant clean separation of energy 
scales enables extraction of the universal scaling behavior 
of dynamics in terms of the gap scale itself; a successful 
description of which is in addition a necessary prerequisite 
for a theory of dynamics and associated transport proper- 
ties at finite-T, which will be considered in a subsequent 
paper. A brief non-technical summary and some conclud- 
ing comments are given in §6. 

2 Background 

The Hamiltonian for the PAM is given in standard nota- 
tion by 

h = E c U> + E( £ / + TflJi-Jflfu, 

+ v£(/£c i(y +h.c.) (2.1) 

i.er 

where the first term describes the uncorrelated conduc- 
tion (c) band with nearest neighbour hopping tij = t. 
The second refers to the f- levels with site energies e / and 
on-site repulsion U, while the final term describes c/f- 
level hybridization via the local matrix element V. We 
focus here on local single-particle dynamics embodied in 
- G{^it) = -i< f(f ia (t)fl) > (and likewise 
G^. CT (u;) for the c-electrons) ; and hence the local spectra 
£>7 i;cr (w) = -7r _1 sgn(w)ImG7 i . (T (a)) (7 = c or f). The key 
feature of DMFT [3-6] is that the (f-electron) self-energy is 
site-diagonal, Sf-.^uj) = SijUf^^); and from straightfor- 
ward application of Feenberg's renormalized perturbation 
theory [36,37], the Gj.^w) are given by 



GUM 



where ui + = uj 



V 2 



V 2 



(2.2a) 

-1 



(2.2b) 



i0+sgn(w) and S 1(T (uj) = S^[{G<; 
the Feenberg (or 'medium') self-energy. Equation (2.2b) 
embodies the connection to a self-consistent impurity model 
that is inherent to DMFT [3-6], since it may be cast in 



the 'single-impurity' form G^. (w) 



Z\ c ff (o>)] _1 with an effective hybridization Z\ c ff(w) = V /2 [w + 
S'i cr (w)] _1 that is to be self-consistently determined. 

In this paper we consider explicitly the symmetric PAM 
with e f = -f , for which n f = JZAfLfia) = 1 and 

n c = J2a( c ia c ia) = 1 f° r au U ■ I n contrast e.g. to the Hub- 
bard model, the problem is thus characterized by two inde- 
pendent dimcnsionlcss parameters: U/t* and V/t*, where 
the hopping is scaled as t = t*/2\/~Z~ c (with Z c — > 00 the 
co-ordination number). 



Victoria E. Smith et al: Kondo Insulators: a local moment approach 



3 



Equations (2.2) are general: independent of lattice type 
and whether/not magnetic ordering arises. Here we con- 
sider primarily the Bethe lattice (BL) (because its spec- 
trum is bounded), for which the Feenberg self-energy 



(2.3) 



We shall moreover focus on the homogeneous paramag- 
netic phase (the antiferromagnetically ordered state is read- 
ily handled but less interesting). In this case equations (2.2) 
reduce to 



G f (o 



LU+- 



■LJ 



V 2 1 

lu+ - E(lo) 4 v ' 



E(co) 



V 2 



co+ - \tlG°(w) 



(2.4a) 



(2.4b) 



written explicitly for the BL; with the conventional single 
self-energy E(lo) = e f + Z f (u>) (= U R (cj)-isgn(cj)U l (uj)) 
defined to exclude the trivial Hartree contribution (of |n/, 
which precisely cancels the bare e/ = — §)■ And the prob- 
lem is particle-hole symmetric, reflected in 

E(lo) = -E(-lo) G^{lo) = -G^(-lu) (2.5) 

with lo = the Fermi level. 

The trivial limits of the model, used below, are two- 
fold. First, V = (for any U), where we denote the 
c-electron Green function by <7o(w), with spectral den- 
sity po(w): for the BL po(ui) = ^-[1 — (w/t*) 2 ] 5 is a 
semi-ellipse of halfwidth t* (from equation (2.4a)), while 
Po{<jS) = [v / 7ri*]~ 1 exp(— [uj/t*] 2 ) for the hypercubic lat- 
tice (HCL) [3-6]. Second, the non- interacting limit U = 
0, denoting the Green functions by g2(oj; V 2 ) with the 
V- dependence explicit (go(ui) = g^(uj;V 2 = 0)). The 
corresponding spectra are related generally (from equa- 
tions (2.2)) by 



d f (u,;V 2 ) = ^d c (u,;V 2 ) 



(2.6) 



and d%(u}:V 2 ) = po(u) — V 2 /lo). The obvious point here 
is that for all V ^ the system is a hybridization gap 
insulator [38]. The gap is soft for the HCL, but hard for 
the BL where the (half) band-gap A q (V 2 ) is given by 



2A° g (V 2 ) = ^tl+W 2 -U 
with corresponding spectrum 



Ud c (u;V 2 ) 



cot* 



(2.7) 



(2.8) 



for A° g < 



\lo\ < A° g + t*. Note trivially that the U = 
spectra cannot be expressed in a one-parameter scaling 
form by suitable dimcnsionless rescaling of lo: no matter 
how lo is thus rescaled (e.g. as w/t*), the g?q each remain 
dependent on the ratio V/t* of bare parameters. 



On increasing U from zero the system remains insulat- 
ing and is perturbatively connected to the non-interacting 
limit, being as such a Fermi liquid (which is wholly com- 
patible with the insulating nature of the state). The limit- 
ing low-w behaviour of the single-particle Green functions 
amounts to a renormalization of the non-interacting limit, 
which is the origin of the renormalized band picture [2,5, 
39]. This follows simply by employing the leading \qw-lo 
expansion of the self-energy S(u>) [2,5], viz 



E(u>) 



1 

Z 



1 



LO 



(2.9) 



from equation (2.5), with Z = [1 - (dE R (uo) / 'Olo)^}' 1 
the quasiparticlc weight (and with E 1 (lo) neglected on the 
grounds that it vanishes in the gap). The leading \ow-lo 
behaviour of the G 7 (lo) then follows (from equations (2.2)) 



as: 



G c (lo)~ 9 c (lo;ZV 2 ) 
G f (Lo)~Zg f Q (Lo;ZV 2 ) 



(2.10a) 
(2.10b) 



Equations (2.10) embody the quasiparticlc behaviour 
of the PAM, akin to the local Fermi liquid quasiparticle 
form for the impurity Green function of the Anderson im- 
purity model (AIM) [2]. While well known per se [2], they 
have an important implication for the scaling behaviour 
of the single-particle spectra -D 7 (u') in the strong cou- 
pling/Kondo lattice regime of large-f/, where the quasi- 
particle weight Z becomes exponentially small [8] (as con- 
sidered in § 5). The renormalized indirect band-gap A g = 
A g (ZV 2 ) is given from equation (2.7) as Z — > by A g = 
ZV 2 /t*, and the spectra D j (lo) follow from equations (2.10, 
2.8, 2.6). Their scaling behaviour in strong coupling fol- 
lows by considering finite lo' = lo/A 9 in the formal limit 
A g oc Z — > 0; and the resultant asymptotic behaviour of 
the Z? 7 (u>) is given by 



UD c (lo) 



1 - 



,12 



V 2 f 4 



1 lo' 2 



(2.11a) 
(2.11b) 



Equations (2.11) show that both t„D c (io) and ^-D f (uj) 
(and not therefore t^D^ {lo)) exhibit one-parameter uni- 
versal scaling in terms of lo' = co/A g , with no explicit de- 
pendence on the bare material parameters U/t* and V/t* 
(which behaviour is naturally not specific to the BL). This 
simple argument does not of course determine either the 
dependence of A g on the bare parameters, or more im- 
portantly the lo'- range in which equations (2.10, 2.11) 
hold: for these a 'real' theory is required. But the ob- 
servation is important in that it suggests the _D 7 (cj) in 
strong coupling should more generally exhibit such scal- 
ing, and provides explicitly the limiting behaviour that, as 
\to'\ = \ui\/A g — ► 0, must of necessity be recovered by any 
microscopic theory; which touchstone will be compared to 
the LMA results in section 5.2 . 
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3 Local Moment Approach (LMA) 

The usual route to single-particle dynamics is via the con- 
ventional single self-energy S(u>). But a determination of 
G(u>) in this way is neither mandatory nor a priori de- 
sirable: in practice, theoretical approaches of this ilk are 
liable to suffer from being essentially perturbative in the 
interaction U (even if self-consistent), thus preventing ac- 
cess to the strongly correlated regime of primary interest. 
The LMA thus avoids such an approach completely and, 
as for the pure impurity models considered hitherto [25- 
33], has three essential elements, (i) Local moments ('/t'), 
regarded as the first effect of interactions, are introduced 
explicitly and self-consistently from the outset: the start- 
ing point is thus static mean-field (MF, i.e. unrestricted 
Hartree-Fock) . This contains two degenerate, local sym- 
metry broken MF states, denoted by a = A or B and 
corresponding respectively to local moments fi = or 

— Severely limited by itself, simple MF nonetheless 
provides a starting point for a non-perturbative many- 
body approach, (ii) The LMA achieves this by employing a 
two-self-energy description that follows naturally from the 
underlying two local saddle points; introducing non-trivial 
dynamics into the two self-energies via their functional de- 
pendence on the broken symmetry MF propagators, (iii) 
The final, key idea behind the LMA is that of symmetry 
restoration: self-consistent restoration of the broken sym- 
metry inherent at pure MF level, and recovery of Fermi 
liquid/quasiparticle behaviour, as discussed below. 

As for the paramagnetic phase of the Hubbard model 
in infinite-d [40], the essence of the approach to the PAM, 
whether at MF level or beyond, is statistical: any given 
site is with equal probability of a =A (ji = |/t|) or B 
(fj, = — type. First consider briefly pure MF, where 
the interaction self-energies are purely static Fock ('bub- 
ble diagram') contributions; given explicitly by = 

— §t/|/i| — ^B-a f° r « = A or B sites. The corresponding 
MF propagators are denoted by g^ai 00 ) (7 = c or 0> an< ^ 
the total Green functions by 



pure MF level the local moment |//| is determined self- 
consistently from the usual MF condition \^\ = \p,(x)\, 
with |/2 1 given by 



|//| 



= f dw [d£ T (w) - d f Al (u>)] (3.3) 

J — 00 



(and dl a (w) = -tt ^gn^Im^w)). 

MF results will be discussed in §4, and their serious 
deficiencies highlighted. We consider now the general case, 
the algebraic structure of which is formally equivalent to 
that at MF level. The full G 7 (w) are expressed as 



where (see equations (2.2, 2.3)) 

V 2 



uj + - S aa (uj) 



— - - A tlG c {io) 



V 2 



(3.4) 



(3.5a) 



(3.5b) 



and the Feenberg self-energy Si a (uj) = jt 2 G c (to) for the 
same reason given above. In contrast to static MF level 
however, the self-energies 



(3.6) 



are now dynamical. They may be separated conveniently 
as Ea<t(u) = —^U\p,\ + E'ao-(w) into (a) the purely static 
Fock contribution — §?/|/t| (that alone is retained at pure 
MF level); together with (b) the dynamical contribution 

Za ct (w) = ^Acr[{5Acr}] that is a functional of the MF 
propagators, and a suitable, naturally approximate choice 
for which (§3.2) determines the extent to which the key 
physics of the problem is captured in practice. The full 
G 7 CT (w) likewise satisfy |/J.-spin symmetry 



(3.1) 



The gZ<j(u) are given for the BL by (see equations (2.2)) 



(3.7a) 



(from equations (3.4, 3.6)); as well as particle-hole sym- 
metry 

Gl a (w) = -Gl_„{-u) (3.7b) 



5a» 



V 2 



uj + + ax 4 



u + + ax 



V 2 



(3.2a) 



(3.2b) 



where x = ^U\fx\; and where (see equation (2.3)) the Feen- 
berg self-energy S^u;) = jt 2 g c (uj) since precisely half the 
(Z c — > oo) nearest neighbours to any given site are of a = 
A (or B) type. The g^^iuj) follow analogously and satisfy 
ffX<rM = 9b-A u ) CT/l-spin symmetry'); while particle- 



hole symmetry is reflected in g^ a (uj) 



-gl-A-u) and 



hence (via equation (3.1)) g 7 (w) = —g 7 (—Lu). And at 



(reflecting S aa {uj) — — £ a _ a (— lo) for the symmetric PAM 
considered). As for the AIMs considered previously [25- 
33], the two-self-energy description inherent to the LMA 
is explicit in equations (3.4, 3.5). The resultant c- and 
f- electron Green functions are, as they must be in the 
absence of an applied magnetic field, rotationally invari- 
ant: Equation (3.4) is correctly a- independent, as follows 
from equation (3.7a). The requisite particle- hole symme- 
try G 7 (uj) = — G 7 (— w) is likewise satisfied, using equa- 
tions (3.4, 3.7b). 

Before proceeding note that, using equation (3.7a), 
equation (3.4) for the G 7 (oj) may be written as G 7 (o>) = 
^ ~^2 a G 7 CT (u;), involving now a spin sum that (from equa- 
tion (3.7a)) is independent of a — A or B (which obviously 
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holds also for the MF Eq. (3.1)). This is entirely equiva- 
lent to equation (3.4). We choose to work with this form 
and, since the a- label is redundant, drop it from now on 
(using a — A implicitly). The basic underlying equations 
on which we subsequently focus are thus 



g 7 M = !E g 2M 



with 



v 2 i 

^G c (co) 

- E a {u>) 4 



E a (io) 



V 2 



and the f-electron self-energies separated as 

S a {ui) = ~U\ji.\+E a {ui) 
with \p\ = \fi(x)\ given by equation (3.3). 



(3.8a) 

(3.8b) 
(3.8c) 

(3.9) 



3.1 E(lu) and Symmetry Restoration 

The conventional single self-energy E(lo) may if desired be 
obtained as a direct byproduct of the two-self-energy de- 
scription inherent to the LMA (although the converse does 
not of course hold). E(lo) is defined by equation (2.4b) for 
G f (u), direct comparison of which to its two-self-energy 
counterpart equations (3.8a, 3.8c) yields 



27M = !(£ T M + AM)- 



[i(AH-AH)] 5 



i(r T H + ^H) 

(3.10) 

where Q{w) = [(G f (w)) _1 + E(lo)}~ 1 is the so-called host/ 
medium f-electron propagator [41]. This relation is general 
(i.e. not specific to the BL), but for the particular case of 
the BL G(io) is given explicitly by 



V 2 



(3.11) 



(using equation (2.4b)). Notice also that E(lo) can equiv- 
alently be defined by equation (2.4a) for G c (lo): direct 
comparison of which to its two-self-energy counterpart 
equations (3.8a, 3.8b) again yields equation (3.10), as it 
must. Given the {E a (w)} and hence (from equations (3.8a, 
3.8b)) G c (co), equation (3.10) enables E(w) to be deter- 
mined. 

The final, important notion underlying the LMA is 
symmetry restoration (SR) [25-33] : self-consistent restora- 
tion of the broken symmetry endemic at pure MF level, 
and correct recovery of the low-u> quasiparticle behaviour 
equations (2.10) that reflects adiabatic continuity to the 
non- interacting limit. As for the AIMs considered hith- 
erto [25-33], whether symmetric or asymmetric, this is 
embodied in the SR condition E^{w = 0) = £i(u> = 0) 



at the Fermi level; and hence E a (uj = 0) = (for ei- 
ther a) for the present particle-hole symmetric problem 
(E a (cj) = i.e. 

Ef{u) = 0) = E^w = 0) - ^U\fi\ = . (3.12) 

If SR equation (3.12) is satisfied, then the leading uj — > 
behaviour of the S a (u>) follows from particle-hole sym- 



metry as (ReE a (u>) =)E^(u 



-{Z a 1 — 1)oj, where 



1 is thus defined and is in- 
= E l a (w) may be neglected 



z a = [i- (dsfiuydu)^] 

dependent of a. (The E\{lo) 
inside a hard gap or, if the gap is soft, because they van- 
ish sufficiently rapidly as w — > 0; which behaviour is guar- 
anteed from the diagrams for E a (u>), §3.2ff.) Using this 
asymptotic behaviour in the basic two-self-energy equa- 
tion (3.8) for the G 7 (w) shows that the quasiparticle form 
equation (2.10) is correctly recovered, with quasiparti- 
cle weight Z = [1 - {dS^i^/duo)^]- 1 = Z a - i.e. the 
leading low-w behaviour of S^- (ui) and U K (tu) coincide, 
£^(u>) = U r (uj) = —(Z^ 1 — l)u> (as may also be veri- 
fied directly from equation (3.10)). And we add that the 
persistence of the insulating gap with increasing interac- 
tion U, implied by the quasiparticle form equations (2.10) 
with Z > 0, is readily shown to be guaranteed only if SR 
equation (3.12) is satisfied. Finally, we note that the SR 
condition (u> = 0) = E± (u> = 0) can be shown to apply 
generically to the paramagnetic phase of the asymmetric 
PAM, which is in general metallic. In this case, as for the 
asymmetric AIM discussed in [30], satisfaction of SR en- 
sures that E 1 {lo) oc lo 2 as ui — > and hence the recovery 
of metallic Fermi liquid behaviour as required by continu- 
ity to the non-interacting limit. Equation (3.12) is thus a 
particular case of the general SR condition, applicable to 
the particle-hole symmetric PAM where the ground state 
is insulating. As for the AIM its imposition — as a single 
condition at the Fermi level lo — — underlies the LMA, 
amounting in practice to a self-consistent determination 
of the local moment \/j,\ (that supplants the pure MF con- 
dition |^| = |/2(x)|, see equation (3.3)); and, most im- 
portantly, generating the low-energy spin-flip/Kondo scale 
that is symptomatic of the Kondo insulators and detailed 
in section 5. 



3.2 Self-energies 

The above discussion is general, and to proceed in prac- 
tice requires specification of the dynamical E a (u>Ys (equa- 
tion (3.9)). The specific class of diagrams retained here is 
naturally motivated on physical grounds. They embody 
self-consistent, dynamical coupling of single-particle exci- 
tations to low-energy transverse spin fluctuations on all 
sites, that is essential in particular to capture the strong 
coupling Kondo lattice regime. The diagrams are shown 
in Figure 1, single lines denoting the MF f-electron prop- 
agators for site i, and wavy lines the local interaction U ; 
the double line propagator denotes the broken symmetry 
host/medium f-electron propagator Qu{lo) specified below 
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(a) { 



(b) 



(c) 




(rr=) 



Fig. 1. (a) Class of diagrams for the f-electron self-energies 
E a (u)) here retained in practice. The interaction U is denoted 
by a wavy line, MF propagators by a single line, and the 
renormalized host/medium propagator (see text) by a dou- 
ble line, (b) Equivalent recasting (including incoming/outgoing 
propagators) to illustrate the spin-flip scattering involved, (c) 
Particle-hole ladder sum in the transverse spin channel; for 
II h (io), spins are reversed. 



(see also [40]). The diagrams may be regrouped as shown 
in Figure 1, which translates to: 



du>i a 
2m 



r^-wijn-^^) (3.13) 



In physical terms these diagrams describe dynamical, 
correlated spin-flip scattering processes: in which having, 
say, added a cr-spin electron to a —cr-spin occupied f-level 
on site i, the —cr-spin hops off the f-level and thus gen- 
erates an on-site spin-flip (reflected in the transverse spin 
polarization propagator 77~ <TCT (u>)); the —cr-spin electron 
then propagates through the lattice/host in a correlated 
fashion, interacting fully with f-electrons on sites j / i 
(reflected in the host/medium Q_ a (ui))\ before returning 
to site i at a later time whereupon the originally added cr- 
spin is removed (which process simultaneously restores the 
spin-flip on site i). The renormalized f-electron medium 
propagator Q- a , which embodies correlated propagation 
of the —cr-spin electron through the lattice, is given ex- 
plicitly for the BL by (cf. equation (3.11)) 



V 2 



. (3.14) 



Its diagrammatic expansion in terms of MF propagators 
and dynamical self-energy insertions E- a (u), and hence 



the infinite set of diagrams implicitly summed in equa- 
tion (3.13) for S a (w), is readily shown to have precisely 
the same topology as its counterpart for the Hubbard 
model employed previously in [40] . We thus refer the reader 
to [40] (section 3.1 therein) for detailed discussion of the 
diagrammatics, here emphasizing just two points, (i) Renor- 
malization of JJ a (u>) in terms of the medium propagator 
Q_ a (u)), rather than G_ a itself, is embodied in the fact [40] 

that G- a , contains dynamic self-energy insertions on any 
site j excluding the original site i (as implicit in the static 
MF interaction on site i appearing in equation (3.14)). 
This accounts in effect for the hard core boson nature of 
the on-site spin-flip, which would be violated if renormal- 
ization in terms of G^_ a was employed (itself necessitating 
inclusion of additional classes of cancelling diagrams), (ii) 
Most importantly we emphasize the strongly renormal- 
ized, and hence self-consistent, nature of the diagrams re- 
tained in E a {u>). This is reflected directly in the fact that 
G- a (uj) (equation (3.14)) depends upon G c (uj) (via the 
'effective hybridization' A eS {w) = V 2 [oj+ - \tlG c {w)]- 1 ), 
which is to be self-consistently determined via solution of 
equations (3.8). 

The transverse spin polarization propagator entering 
equation (3.13) for E a {uj) is given at the simplest level, 
shown explicitly in Figure 1, by an RPA-like particle- hole 
ladder sum in the transverse spin channel; viz 



n a ^(cu) = °n a ^(uj) [i - u°n a ^(Lu)]~ 



(3.15) 



where the bare particle-hole bubble is itself expressed in 
terms of the broken symmetry MF f-electron propaga- 
tors {gl}- Our subsequent discussion refers explicitly to 
n a ^ T (co) described at this level; in §5 a further renormal- 
ization is also considered, wherein °77 <t ^ t (cj) and hence 
II cr ~ cr (Ld) are expressed in terms of the fully self-consistent 
host/medium propagators {Ga} (i-e. with all propagators 
in Figure 1 renormalized in terms of the double line Q a 's) . 
Whichever level is employed, the a n a ~ cr and hence 
are readily shown to be related by n~ c " T {w) = 7T CT ^ T (— u>); 
whence only one such, say 77 ^ (u>), need be considered ex- 
plicitly. Using this, and the Hilbert transform for 77 H (u), 
equation (3.13) reduces to the following form convenient 
for later analysis (§5.1) 



£ T (w) = U 2 



duji 
ir 



Im J7+- ( Wl ) [6»(wi )g^ (wi + w) 
+ (3.16) 



(with Ei(ll>) = — w)); where 9(x) is the unit step 
function and Ga = Qt + ^<r i s separated into the one- 
sided Hilbert transforms 



Q${w) = / dwi 



uo — uii ± i0 + 



(3.17) 



with D a (uj) = — 7r _1 sgn(a;)Im(7 -(ti)) the corresponding spec- 
tral density. Note also that E a (uj) depends explicitly on 
U (via the interaction vertices), as well as implicitly on 
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x = \U\ii\ via the dependence of the MF propagators 
thereon (equations (3.1, 3.2)). In particular, the symmetry 
restoration condition equation (3.12) has the functional 
form 

Zf(LJ = 0;U,x) = ±U\ji(x)\ (3.18) 
exploited below. 

3.3 Practice 

We now draw together the preceeding discussion and spec- 
ify what we find to be an efficient practical algorithm to 
solve the LMA-DMFT equations. This consists of the fol- 
lowing basic steps: 

(i) 'Startup'. For any given x — ^U\/i\, equations (3.1, 
3.2) are first solved to give the MF propagators g 7 (w) 

(= IL^H with s2M = slaM and t = c or f )- 

From this the polarization propagators n a ~ a (u)) follow 
(equation (3.15)); and with G c (uj) ~ <? c (w), the initial f- 
electron medium propagator Q_ a (u)) follows on comparing 

equations (3.14, 3.2b) as G-tr{w) — 9- CT ( w )- The startup 
E a (uj) = E a {uj; U, x) follows directly from equation (3.13). 

(ii) The symmetry restoration condition S^ {lo = 0; U, x) = 

7^U\p,(x)\, is then solved for U given x = ^U\fj,\ (or vice- 
versa which, while entirely equivalent, is less efficient in 
practice); the local moment \/j,\ follows immediately. 

(iii) With the resultant = —^U\p,\ + S a {io) (equa- 
tion (3.9)), equations (3.8a, 3.8b) are then solved directly 
for G c {uu) (with G?(uj) following in consequence from equa- 
tions (3.8a, 3.8c)). 

(iv) The resultant G c {uj) is then used in equation (3.14) 
to obtain a new and hence via equation (3.13) 
a new S a {un). Now return to step (ii) and iterate to self- 
consistency. 

We find the above algorithm to be efficient, typically 
converging after 4-5 iterations and computationally fast 
on a PC. It is easily generalized to the case (§5) where the 
77 ^ 's are renormalized in terms of the {^ CT }; as well as to 
encompass the hypercubic lattice, section 5.2 (or indeed an 
arbitrary lattice DOS po( w ) (see §2) that may be employed 
in materials modelling applications of DMFT). 

Results from pure MF (step(i) alone, with \fj,\ deter- 
mined from \n\ = \ jl(x)\) will be discussed in §4. We shall 
also discuss separately (§5.1) '1-loop' results obtained from 
a single iterative loop of the above scheme (steps (i)-(iii)): 
this has the advantage of being analytically tractable, and 
contains much of the key behaviour found in the full iter- 
ative solution (§5.2). 

Before proceeding we mention the issue of stability 
that is important for the PAM, as for the AIMs [25-33] 
within the LMA: the fact that Rc77 + ~(w = 0) > of 
necessity (as follows directly from its Hilbert transform). 
For this to be satisfied, using equation (3.15) (and that 
lm°n+-(uj = 0) = 0), < URe°n+-(cu = 0) < 1 is re- 
quired. An explicit expression for Re /! -1 (0) is however 
readily obtained (see e.g. [25,30]), viz 

(7Rc°7T+-(a> = 0) = ^M (3.19) 



with \j2{x)\ given by equation (3.3) (and x = ^U\/j,\). For 
stability, \n\ > \fi(x) \ is thus required. The broken symme- 
try pure MF solutions, for which = \p.(x)\ determines 
the local moment, are thus properly stable (always, which 
we note would not be the case if restricted Hartree-Fock, 
with |^|=0 enforced a priori, was employed). But they lie 
on the 'stability border', with £/Re°Z7+-(w = 0) = 1. The 
latter in turn implies, from equation (3.15), that the trans- 
verse spin propagator 77 ^ (u>) contains a pole at u = 0; 
reflecting physically the fact that the pure MF state is, 
locally, a degenerate doublet. The latter behaviour is cor- 
rect for a local moment phase, which for the PAM means 
the zero-hybridization limit where the f-electrons decou- 
ple from the conduction band (which limit we add the 
LMA recovers exactly). It is not however correct for the 
Kondo insulating state that is adiabatically connected to 
the non-interacting (singlet) limit: here by contrast the 
characteristic energy scale for the local spin-flips is non- 
zero, and on the order of the Kondo scale that typifies the 
Kondo insulators. 

The above behaviour is however entirely specific to 
the pure MF level of self-consistency, i.e. arises only if 
\n\ is determined by \/j,\ = \fi(x)\. And the central point 
is that within the LMA the local moment is determined 
from the symmetry restoration condition equation (3.18) 
(step (ii) above in which \/j,\ > \fi(x)\ is always found). In 
consequence Im77 + ~(u>) contains not an ui = spin-flip 
pole, but rather a resonance centred on a non-zero fre- 
quency uj m . This is the low-energy scale characteristic of 
the Kondo lattice, its origin within the LMA thus stem- 
ming from self-consistent imposition of symmetry restora- 
tion; and its physical significance being that it sets the 
timescale r ~ h/u) m for restoration of the locally broken 
symmetry /degeneracy inherent at pure MF level. We add 
moreover that for V ^ symmetry restoration is found to 
be satisfied for all U > 0. Its breakdown at finite U would 
signal an underlying quantum phase transition (such as 
arises in the pseudogap impurity model, see e.g. [31-33]); 
which, in concurrence with general belief, does not there- 
fore arise for the paramagnetic phase of the PAM within 
DMFT. 



4 Mean-field 

Our discussion of pure MF is brief: it is fundamentally 
flawed, lacking even an insulating gap beyond small inter- 
action strengths, as illustrated below. Understanding its 
deficiencies is however important, for it ultimately under- 
pins a successful description of the Kondo insulating state. 
Moreover as demonstrated in §5, the primary deficiencies 
of pure MF arise on energy scales on the order of the insu- 
lating gap or multiples thereof. Such scales are naturally of 
paramount importance, but arc nonetheless small — ex- 
ponentially so in the strong coupling Kondo lattice regime 
— compared to bare ('non-universal') scales on the order 
of U, i* or V. This suggests that MF might provide an es- 
sentially sound description of dynamics on non-universal 
energy scales; an issue considered in §5. 
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co/t t co/t, 



Fig. 2. Illustrative mean-field results for V 2 ft* = 0.2. (a) MF local moment \fio\ vs. U/t*. MF f-electron spectra, t*dJ (ui) vs. 
uj/U for (b) x = \U\no\ = 0; (c) x = 0.9 (U/U ~ 2.1); (d) a; = 2 (U/U ~ 4.2). 



Figure 2 shows the t//i* dependence of the MF mo- 
ment for V 2 /t 2 = 0.2, determined as in §3 from = 
\p.(x)\ (x = ^U\n\) and denoted by |/xo ! - The local mo- 
ment first becomes non-zero for U/U ~ 1.41, increasing 
rapidly thereafter with U/t* such that for U/t* > 2 or so 
the moment is well formed and close to saturation. Such 
behaviour sets in at even lower U/t* with decreasing hy- 
bridization V 2 /t 2 , reflecting the atomic limit incipient as 
V — > where the f-levels decouple from the band. 

Illustrative spectral evolution with U is given in Fig- 
ures 2(b)-(d), where the f-level spectrum t*d^(u)) vs. uj/U 
is shown for x = \U\p Q \ = 0,0.9 (U/U ^ 2.1) and 2 (U/U ~ 
4.2). The characteristics of the non-interacting hybridiza- 
tion gap insulator (Figure 2(b)) are rapidly lost with in- 
creasing U: MF produces an insulator-metal transition 
occurring at U/U ~ 1.48 just after local moments form 
(and arising generally at x — 2V 2 /t„ as follows from equa- 
tions (3.1, 3.2a) at the Fermi level u> = 0). With increasing 
U/U, Figure 2(c, d), Hubbard satellites at |w| ~ f form 
in the f-electron spectra, rapidly acquiring dominant spec- 
tral weight at the expense of intensity on the 'band' scales 
|w| < i*; the converse occurring for the MF c-electron 
spectrum d c (uo) (not shown) where spectral intensity on 
the band scales increases with U/t*. 

The simple nature of the MF spectra in strong coupling 
(U 3> 2V 2 /t*) is in fact seen directly from equations (3.1, 
3.2). For x = lU\m,\ ~ t > V 2 /U, and for \w\ < 
equations (3.1, 3.2a) yield 



d c (ui) ~ Po{<-o) characteristic of the decoupled V = limit 
(a semi-ellipse for the BL considered explicitly). In prac- 
tice this behaviour is well attained for x > 2 or so, and 
in consequence the f-electron g£(u>) follows from equa- 
tion (3.2b) as gf(u) ~ [lo+ + erf - V 2 g (uj)] ~\ The re- 
sultant dl(oj), while dominated by the satellites at \u>\ ~ 
f , nonetheless contains a low-intensity continuum on the 
band scales \u>\ < U, given explicitly by: 



u 

2 



oj + \tlg c (u) 



= go(uj) (4.1) 



(with go(u>) the V" = conduction band Green function, 
§ 2); i.e. V and U drop out, producing a conduction band 



V 2 



(4.2) 



This is evident in Figure 2; it will also prove important in 
determining (§5.1) the Independence of the Kondo scale 
within the LMA. 

The 'insulator-metal transition' is of course entirely an 
artifact of MF, and the resultant 'metal' not a Fermi liq- 
uid: the pure MF single self-energy Zmf(w) = ^mf^) — 
isgn(w)Z , ^ [F (cj) is obtained from the general result equa- 
tion (3.10) using S a = —ax; is given by 



(4.3) 



(which we note is not purely static), and is such that 
E l MF (ui = 0) ^ at the Fermi level if d c (ui = 0) oc 
Img c (aj = 0) is non-zero — thus violating Fermi liquid be- 
haviour. But the perversity of MF in predicting a metallic 
state for the symmetric PAM which describes the Kondo 
insulator, does not survive the inclusion of correlated elec- 
tron dynamics in the self-energies S a (uu). To which we now 
turn. 
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Fig. 3. LMA spectra, all-scales overview. Top panels: conduction electron spectra; bottom panels: f-electron spectra, (a) For 
x — ^U\fi\ = 1 (U = 2.2), (b) For x = 3 (U = 6.1); with V 2 = 0.2 in either case. Corresponding results at pure MF level are 
shown for comparison (dashed lines). 



5 Results 

We begin with a brief overview of spectral evolution on 
all energy scales, and from weak to strong coupling in- 
teraction strengths, before turning (§5. Iff) to the central 
behaviour at low-energies, and in the strongly correlated 
regime. 

With t* = 1 as the unit of energy, and V 2 = 0.2, 
Figure 3 shows representative c- and f-electron spectra 
obtained from the LMA: for x — |f7|/i| = 1 an d 3 (corre- 
sponding respectively to U = 2.2 and 6.1). Corresponding 
results at pure MF level are shown for comparison. The 
x = 1 example is transitional between weak and strong 
coupling behaviour. Here the Kondo insulating gap, which 
the LMA preserves correctly for all U, is discernible (just) 
on the scales shown; and Hubbard satellites, apparent in 
both Df(u>) and D c (oj), are just 'breaking away' from the 
main band. The x = 3 example by contrast is typical of 
strong coupling behaviour (which in practice is reached by 
x ~ 2): the gap is exponentially small (§5.1, 5.2) and as 
such not resolved on the scales shown; and the satellites 
are well formed, dominating the net spectral intensity in 
D* (oj) (as is physically obvious) and giving an ever dimin- 
ishing contribution to D c (lo). 

Comparison to the MF spectra is revealing, particu- 
larly for the strong coupling example x — 3 (Figure 3(b)). 
The Hubbard satellites are, unsurprisingly, broadened and 
shifted slightly from their MF counterparts. However for 
the c-electron spectrum it is clear that on the band scales 
\u)\ < t», and excepting the (all important) low-energy 
gap region, the LMA and MF spectra are nigh on coinci- 
dent; supporting the notion that MF itself provides a rea- 
sonable description of the conduction band on bare/non- 



universal energy scales. We emphasize that this applies 
only to D c (ui) and not to the f-electron spectra: here, as 
seen directly from Figure 3, the MF spectra are deficient 
throughout the band, reflecting the dominance of corre- 
lated electron dynamics for the f-levels, that are simply 
absent in MF. 

The LMA results shown in Figure 3 refer to the fully 
self-consistent level described in §3.3. We add however 
that on the full energy regime shown these differ insignifi- 
cantly from the LMA results obtained either at the simpler 
1-loop level (§3.3), or with the transverse spin polariza- 
tion propagators 7T CT ~ CT further renormalized in terms of 
the host/medium propagators {^o-}; as illustrated in Fig- 
ure 4 for x = 3 (U = 6.1). The differences between these 
different levels arise of course on the low-energy gap scale, 
as pursued in the following sections. 

The comments above refer to intermediate to strong 
coupling behaviour, with a natural emphasis on the lat- 
ter. But what of weak coupling? Here, as for the AIMs 
considered hitherto [25,30,31], the LMA is readily shown 
to be perturbatively exact to/including second order in 
the interaction U. Figure 5 demonstrates the point, show- 
ing LMA results for D f (u) with V 2 = 0.2 and U = 1 and 
0.25, compared to those from simple second order pertur- 
bation theory in U (SOPT): with decreasing U the LMA 
spectrum clearly reduces to that from SOPT, and for the 
lower U shown the two are essentially indistinguishable on 
all energy scales. 

We turn now to the central issues: low-energy dynam- 
ics on the scale of the Kondo insulating gap, particularly 
in the strong coupling regime; and, most importantly, uni- 
versal spectral scaling in terms of the gap itself. 
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Fig. 4. All-scales comparison of LMA spectra at different levels 
of description, for x = ^U\(jl\ = 3 (U = 6.1) with V 2 = 0.2. Top 
panel, c-electron spectrum; bottom, f-electron spectrum. Solid 
line: self-consistent LMA ('LMA I', see text, §5.2). Dashed line: 
with further renormalization of IJ"^ 7 propagators as described 
in text ('LMA IP). Dotted line: 1-loop LMA results. For the 
f-electron spectra, the three are barely distinguishable on the 
scales shown. 




Fig. 5. Weak coupling LMA f-electron spectra (solid lines), for 
U — 1 and 0.25 (with V 2 = 0.2). Corresponding results from 
second order perturbation theory in U (SOPT, dashed lines) 
are also shown; for the lower-?/ they are indistinguishable from 
the LMA. 

5.1 1-Loop LMA 

We begin by discussing the LMA at 1-loop level, viz a 
single iterative loop (steps (i)-(iii)) of the general scheme 
discussed in §3.3. Results from the full iterative solutions 
will be considered in §5.2; but the simpler 1-loop level is 
itself important because it offers an analytical handle on, 
and insight into, the problem. 

The inset to Figure 6 shows the low-energy behaviour 
of the 1-loop f-electron spectrum vs. lo(= w/i*), i.e. on an 
'absolute' scale; for V 2 = 0.2 and two different interaction 



strengths, U = 6.1 (x = \U\pi\ = 3) and U = 8.1 (a; = 4). 
The spectra are quite distinct and clearly dependent on 
the bare material parameters, with the gap narrowing 
strongly under the modest increase in U. The main part of 
Figure 6 by contrast shows 2ttV 2 D^ (lo) vs. lo' = lo/[ZV 2 ] 
where Z = [l-(dE*(w) / '0w)u,=o] _1 (§3.2) is the quasipar- 
ticle weight. The point is obvious: the spectra collapse to 
a common form; which universal scaling is indeed seen to 
be in terms of the gap scale A g = ZV 2 , as required from 
the quasiparticle form discussed in §2 that embodies per- 
turbative continuity to the non-interacting limit (and full 
comparison to which is given in §5.2, see Figure 9 below). 
The scaling spectrum V 2 D^(lo) is entirely independent of 
the two bare parameters U and V , which enter solely via 
the dependence of the gap scale thereon; and we empha- 
size that this applies to V 2 D^(lo) and not therefore to 
Df(u))[= t^D^ (lo)] (i.e. if the calculations of Figure 6 are 
repeated for different V 2 , (lo) vs. lo' is not itself uni- 
versal but V 2 D i is). Before proceeding we add that, as 
for the pure impurity models considered hitherto [25-33], 
the small spectral 'dip' at \lo'\ ~ 2.5 is entirely an artifact 
of the specific RPA-like form for the n a ^ J (io) employed 
here; it can be eliminated [26] but we are content to live 
with it in the following. 

The first obvious question is how the quasiparticle 
weight, and hence A g = ZV 2 , depends upon the bare pa- 
rameters in the strong coupling regime of interest. This 
can be determined analytically as now outlined, focus- 
ing on equation (3.16) for S^(lo). At the 1-loop level the 

medium propagator Q± in equation (3.16) reduces to the 
MF propagator as explained in section 3.3, step (i). 
But in strong coupling the spectral weight of Imi7 H (lo) 
is confined to lo > 0, with J °°(cL;/7r)Imi7 H (lo) = 1; be- 
haviour that reflects physically the strong coupling satu- 
ration of the local moment (|/z| — > 1). As discussed in §3.3 
the resonance in Imi7 H (lo) is (by definition of Lo m ) cen- 
tred on the low-energy spin- flip scale Lo m ; and on scales of 
this order ~ (lo) is slowly varying, whence equation (3.16) 
for U^(lo) reduces asymptotically to 

Sf(Lo) ~ U 2 Reg{-(LO + Lo m ) (5.1) 

But Rcg^~ (u>) is given by the one-sided Hilbert transform 
(equation (3.17)); which as lo — > is dominated by the log 
singularity arising because (§4) the MF spectrum d£(w) in 
strong coupling is itself non-vanishing at the Fermi level 
lo = 0. This lo — > behaviour is captured by 

Re fl f-(w) ~ d{(0) y° 

= 4(0)lnQ£) (5.2) 

where a high energy cutoff of order t* is employed (its 
precise value being immaterial, the important point be- 
ing that the prcfactor to the log divergence is precisely 
<f[(0)). But in strong coupling ^(lo = 0) is given from 
equation (4.2); using which equations (5.1, 5.2) yield the 
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Fig. 6. 1-Loop LMA results for scaling behaviour of f-electron spectrum. Inset: 2tt^—D^{lu) on an 'absolute' scale vs. w/i«, 
for U/t* = 6.1 (solid line) and 8.1 (dotted line). Main figure: collapse to scaling form (solid line), same results shown vs. 
iot t ./(ZV 2 ) = uj/Ag. The scaling spectrum for the Anderson impurity model (see text) is also shown in the main figure (dashed 
line) . 



desired low-frequency behaviour 
Ef(u) 4F 2 p (0)ln 



UJ. n 



(5.3) 



(with po(0) = 2/(7rt*) for the BL, although we add that 
equation (5.3) holds for an arbitrary underlying lattice). 
From this the quasiparticle weight 



[i-(dsf(u)/du) u 

for io m /V 2 p (0) - by 



follows directly, being given 



V 2 

0* 



(5.4a) 
(5.4b) 



with the latter explicitly for the BL. 

The gap and spin- flip scales are thus equivalent. And 
their mutual dependence on bare material parameters fol- 
lows from symmetry restoration (§3.3, step (ii)), viz from 
equation (3.12) in strong coupling (where \p\ — > 1) via 
= 0) = \\J\ which, combined with equation (5.3) 
gives the desired result 



i*exp 



(— 



U 



\%V 2 Po (0) 



= i*exp 



16V 2 



(5.5a) 



(5.5b) 



showing that A g = ZV 2 /i* is indeed exponentially small 
in strong coupling. 

Once symmetry has been restored as above, the self- 
energies E a (u)) (equation (3.9)) for all ui follow imme- 
diately. The final step in the 1-loop analysis (step (iii), 
§3.3) is to take the resultant E a (uj) and use them in equa- 
tions (3.8a, 3.8b) to solve self-consistently for G c {lo). This 



in turn determines directly the effective f-electron hybridiza- 
tion ZWM = V 2 [io+ - 5(c)]- 1 = V 2 [io+ - i^ G c (w)] -i 

(BL); in terms of which G' (ui) and hence D* (ui) follow 
without further ado from equations (3.8a, 3.8c), viz 



(5.6) 



The 1-loop results shown in Figure 6 have of course been 
obtained in this way, 

There is however a simpler, certainly crude but nonethe- 
less revealing approximation that may be employed at this 
final (step (iii)) stage to determine G j F (cu): namely to re- 
place the G c -dependence of A(oj) = A[G C ] by A(uj) ~ 
A[g c ] in terms of the MF propagator (as used also in 
the 'startup' step (i) (§3.3), there leading to Q- a ~ gL a 
in S a {uj)). With this, the resultant A(uj) — A^{lo) — 
isgn(aj) A\{oj) reduces in strong coupling, using Eq. (4.1), 
toA(u) w V 2 g Q {uj) withff (w)(= Re# (w)-isgn(u;)7rpoM) 
the V = conduction electron propagator (a result read- 
ily shown to hold for a general lattice). In the strong 
coupling regime where A g cx Z — ► 0, the w-dependence 
of go(w) is moreover irrelevant, occurring as it does on 
non-universal scales on the order of t*; whence A{lo) w 
— isgn{uj)iTV 2 po{Q) and equation (5.6) becomes 

G f (cu) « i ^ [w+ - S a (u) + isgn{aj)Ao\ (5.7) 

a 

where A = irV 2 p o (0). 

Equation (5.7) is simply the local impurity Green func- 
tion for a metallic Anderson impurity model with hy- 
bridization strength A = irV 2 p (0) = 2V 2 /t* (BL). The 
resultant impurity scaling spectrum TrZio-Dimp^) is also 
compared to the 1-loop ttAqD^ (oj) in Figure 6. The AIM 



12 



Victoria E. Smith et al: Kondo Insulators: a local moment approach 



spectrum is of course metallic (with 7rZ\ Mmp(^ = 0) = 1 
as required by the Friedel sum rule [2]), in contrast to 
the gapped PAM case; but beyond the gap the spectral 
tails of the two rapidly coincide. This is natural, for the 
nature of the 1-loop calculation is clear: the 1-loop self- 
energies E a {uj) are themselves those of the underlying 
AIM. And since the associated low-energy scale is deter- 
mined from symmetry restoration via S^(uj — 0) = 0, 
the low-energy gap scale is simply the Kondo scale for 
the AIM itself; as seen directly from equation (5.5) ex- 
pressed as A g jt* ~ exp(— nU/8A ), recovering the exact 
exponent for the symmetric AIM [2] . The gap in the PAM 
spectrum (which is strictly soft at 1-loop level, albeit not 
visibly so in Figure 6) arises from self-consistent solution 
of equations (3.8a-c) for G c {uo) and hence G^(uj), given 
the effective AIM S a (uj). In that sense the 1-loop level is 
akin to an average t-matrix approximation [20], but for- 
mulated within the two-self-energy description inherent to 
the LMA and with symmetry restoration ensuring Fermi 
liquid behaviour (in the general sense of adiabatic conti- 
nuity to the non-interacting limit). 

Despite the relative simplicity of the 1-loop LMA, the 
relation of the AIM implicit therein is both natural and 
substantially correct: NRG calculations for example [8] 
(see Figure 1 therein) indeed show, beyond the gap scale, 
a close similarity between the PAM {uu) and the AIM 
spectra. Moreover since the spectral 'tails' of the 1-loop 
PAM and AIM are common (Figure 6), recent LMA re- 
sults for the latter [26] may be used directly to infer their 
analytical form for V 2 D i (w) at 1-loop level; specifically 

w vo(o)^H^{ [(4/7r)ln( ;_ l)]2+1 



[(4/7r)ln(|i>|)] 2 + 25j V • ' 

where Co — uj/uj m (and which form is known [26] to agree 
quantitatively for |a>| > 5 with NRG results [32] for the 
AIM itself). These slow logarithmic tails persist in the 
V 2 Df (lu) scaling spectrum beyond 1-loop level (§5.2) and, 
as for the AIM [26,29], are ultimately responsible for the 
'high' temperature logarithmic behaviour of transport prop- 
erties such as the resistivity. 

In physical terms the limitations of 1-loop level are 
nonetheless self evident; it misses the dynamical inter- 
site correlations embodied (as discussed in § 3.2) in the 
functional dependence of the self-energies E a {uj) upon the 
renormalized f-electron medium propagators {Go} (which 
at 1-loop level are simply replaced by their MF counter- 
parts). 



5.2 Fully self-consistent LMA 

The full LMA solutions are now considered. In the fol- 
lowing, we denote by 'LMA I' results obtained from the 
class of diagrams contributing to E a (oj) shown explicitly 
in Figure 1; while 'LMA IF refers to those obtained from 
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Fig. 7. Dependence on bare parameters of the low-energy scale 
uj m oc A g = ZV 2 /t„, arising from the LMA at 1-loop (circles), 
LMA I (squares) and LMA II (triangles) . 

the further rcnormalization of H a ^ a in terms of the f- 
electron medium propagators (i.e. all propagators in Fig- 
ure 3 renormalized in terms of the double line CJ CT 's). 

Figure 7 shows the dependence on bare parameters of 
the low energy scale u> m oc A g = ZV 2 /t* arising from 
LMA I /II for the Bethe lattice, in comparison to that ob- 
tained at 1-loop level (Eq. (5.5), which form is confirmed 
in Figure 7, with the prefactor <~ 0.2i* determined nu- 
merically). Two points should be noted here. First that 
the lattice scale is enhanced over its 1-loop counterpart 
(which as above is equivalently that of the underlying 
AIM); qualitative behaviour that is also found in NRG [8] 
and QMC [10] calculations for the symmetric PAM. The 
second point is the exponential dependence of the scale, 
u> m oc cxp[— XU/8V 2 po(0)] with A = 1 at 1-loop level. 
From Figure 7 we find that A = 1 remains at the level of 
LMA I, while A ~ 0.85 is found from the fully renormal- 
ized LMA II. The latter is in qualitative agreement with, 
but somewhat larger than, NRG calculations [8] which 
yield A ~ 0.7; which in turn are comparably in excess 
of A = i arising from the Gutzwiller approximation [23]. 
The dependence of the low-energy scale on bare model 
parameters is however a subsidiary issue compared to the 
scaling behaviour of dynamics in terms of the gap scale 
itself. The latter is illustrated in Figure 8 which shows the 
f-electron scaling spectra (u) vs. u>' — u/A g , for 

both LMA I and II; the corresponding c-electron spectra 
t*D c (u) are shown in the left inset. Despite the differences 
in the absolute values of the low-energy scale (embodied 
in Figure 7), the resultant LMA I/II scaling spectra differ 
very little from each other; and exhibit the same qualita- 
tive behaviour as found at 1-loop level (Figure 6), includ- 
ing the logarithmic spectral tails in V 2 D^(uj) evident in 
the right inset to Figure 8. 

As discussed in §2, adiabatic continuity to the non- 
interacting limit requires that on sufficiently low energy 
scales the resultant spectra should conform to the dic- 
tates of Fermi liquid theory, as reflected in the quasi- 
particle form equations (2.10, 2.11). In Figure 9 the f- 
and c-electron scaling spectra arising from LMA II are 
compared directly to the limiting quasiparticle behaviour 
equation (2.11) (corresponding comparisons for LMA I 
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Fig. 8. f-electron scaling spectra 2-K^—D f (u>) vs. u)/A g arising from LMA I (dashed line) and LMA II (solid line). Right inset: 
on an expanded scale to show the logarithmic tails. Left inset: corresponding c-electron spectra from LMA I/II. 



and the 1-loop results are very similar, as evident from 
Figures 6,8). Agreement with the quasiparticle form is es- 
sentially perfect close to the gap edges; and in practice is 
followed quite closely up to — \uj\/A g ~ 3 or so. Be- 
yond this however, the quasiparticle form fails to capture 
the logarithmic tails of the scaling spectrum, decaying in- 
stead as V 2 Df(uu) ~ 1/K| 2 (equation (2.11b)). This is 
of course natural since the quasiparticle behaviour is con- 
fined strictly to the limiting low-a/ behaviour, and it mir- 
rors the situation arising in the metallic AIM where the 
quasiparticle form is a trivial Lorentzian [2] that likewise 
fails to capture the slow spectral tails of the Kondo reso- 
nance [26]. 



Finally, an obvious question arises: how strongly do 
the scaling spectra depend upon the host lattice? For the 
metallic AIM the answer is not at all: in the strong cou- 
pling Kondo regime the w-dependence of the host density 
of states to which the impurity is coupled is entirely irrel- 
evant (see e.g. [2]), and the scaling form of it A D- lrav {uj) 
is host-independent. But this is not of course the case for 
the PAM. The preceeding results have been given explic- 
itly for the BL where the free (V = 0) conduction band 
po(w) is semielliptic. To illustrate the influence of the lat- 
tice, Figure 10 compares the LMA I scaling spectra for 
the BL with those for the hypercubic lattice (HCL, where 
po(w) is an unbounded Gaussian). While the Kondo insu- 
lating gap for the BL is hard, it is strictly soft for the HCL, 
albeit that this is barely visible in Figure 10 since the re- 
sultant Z? 7 (w) are exponentially small close to the Fermi 
level. But beyond the immediate vicinity of the gap, and 
despite the very different nature of the host po( w )'s, the 
scaling spectra are seen to be qualitatively very similar; 
supporting the view that for local dynamical properties 
the one-electron 'details' of the lattice play but a minor 
role. 
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Fig. 9. Comparison of LMA II scaling spectra (solid lines) to 
the limiting low-energy quasiparticle form (dashed lines). Top 
panel, c-electron scaling spectra; bottom, f-electron spectra. 

6 Discussion/Outlook 

We have developed in this paper a local moment approach 
to T = single-particle dynamics of the symmetric PAM, 
the basic microscopic model for understanding Kondo in- 
sulating materials [34,35]. The necessary criteria for a suc- 
cessful description of the problem appear to be met by the 
LMA, handling as it does all energy scales and interaction 
strengths while satisfying the requirements of Fermi liquid 
theory at sufficiently low energies. 
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Fig. 10. LMA I scaling spectra for the hypercubic lattice (solid 
lines), compared to those for the Bethe lattice (dashed). 



Particular attention has been given for obvious phys- 
ical reasons to the strong coupling (i.e. large U) Kondo 
lattice regime, believed to be appropriate to the Kondo 
insulator materials mentioned in §1. We have shown that 
within DMFT+LMA one nevertheless recovers correctly, 
on sufficiently low energy scales, an 'insulating Fermi liq- 
uid' behaviour which evolves continuously from the non- 
interacting hybridization gap insulator. For example, the 
scaled conduction electron and f-electron spectra have a 
renormalized non-interacting form (Eqs. (2.10, 2.11)) for 
\w\/A g < 2 — 3 (see Fig. 9); where A g , the renormalized 
gap, is reduced from the non-interacting hybridization gap 
by the quasiparticle weight factor Z. For larger en- 
ergies the (scaled) spectra deviate rapidly and substan- 
tially from the non-interacting form, but nonetheless re- 
main characterized by a single low-energy scale A g oc Z 
that is exponentially small but in general enhanced over 
its counterpart in the dilute (AIM) limit [42] . The spectra 
thus exhibit 'universal scaling' in terms of co/A g (Fig. 9), 
which we find to be dominated by slow logarithmic tails 
just as for the metallic AIM, and for which the LMA pro- 
vides analytic results (see e.g. Eq. (5.8)). We naturally 
expect that these features will also show up in the dy- 
namic and transport properties of the model at finite-T, 
a subject to which we will turn in a subsequent paper. 

Finally, the present work has been confined intention- 
ally to the symmetric PAM relevant to the Kondo insula- 
tors. This is of course a special, albeit physically impor- 
tant limiting case of the asymmetric PAM with arbitrary 
conduction band filling n c , which for n c ^ 1 describes the 
metallic heavy fermion compounds; and extension of the 
LMA to encompass which is currently in hand. 
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